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Any  one  who  considers  arithmetical  methods  of  pro- 
ducing random  digits  is,  of  course,  iii  a state  of 
sin. 

John  von  Neumann 


He  that  is  without  sin  among  you,  let  him  first 
cast  a stone  ... 

John,  8:7 


I 


I.  Introduction  - Why  so  many  tests 

Because  of  its  applications  in  simulation,  statistical  problems, 
and  some  areas  of  numerical  analysis,  an  indispensable  piece  of  soft- 
ware on  a computer  today  is  a copious,  fast  and  reliable  source  of 
uniformly  distributed  random  numbers.  Early  in  computer  history  the 
idea  of  inputting  a table  of  random  numbers  was  rejected  and  the  in- 
dustry moved  to  a programmed  arithmetic  process  for  producing  sequences 
of  random  numbers.  If  the  definition  of  random  means  loosely  "each 
number  obtained  by  chance,  having  nothing  to  do  with  other  numbers," 
then  quite  obviously  numbers  generated  in  such  a deterministic  way  can 
only  appear  to  be  random  and  may  be  termed  "pseudo- random". 

The  most  popular  and  successful  random  number  generators  today  use 
a linear  congruential  method  which  depends  on  the  word  size  of  the  com- 
puter. A program  RANDUL  exists  on  the  SPL  and  SYS.OLIB  of  the  ASC 

Note:  Menutcript  submitted  April  29,  1977. 
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which  implements  this  widely  used  industry  algorithm.  Knuth,  however, 
warns  of  the  pitfalls  of  blindly  accepting  handed  down  algorithms  for 
random  number  generators^  and  Coveyou  and  MacPherson^  warn  of  insuf- 
ficiencies with  using  this  algorithm  and  a 32-bit  word  size. 

A second  generator,  VRANF,  is  also  available  on  the  SPL  and 
SYS.OLIB  of  the  ASC  which  uses  a variation  of  the  linear  congruential 
algorithm  and  the  computer's  vector  speed  to  generate  its  pseudo-random 
numbers.  It  is  also  a common  fallacy  that  one  can  take  a good  gener- 
ator, modify  it  a little  and  get  an  even  more  random  sequence.® 

These  warnings  raised  the  serious  need  for  thorough  testing  of 
both  algorithms.  The  theory  of  statistics  provides  us  with  quantitative 
measures  of  randomness;  however,  there  is  no  end  to  the  number  of  sta- 
tistical tests  that  can  be  performed.  Many  tests  were  performed  kno^'ing 
that  although  not  sufficient  to  prove  randomness,  every  test  passed 
built  more  and  more  confidence  in  the  random  sequence.  Thus,  the  first 
set  of  questions  the  paper  attempts  to  answer  has  to  do  with  "goodness", 
(a)  Do  the  two  programs  generate  sequences  which  behave  as  if  random  or 
more  simply  are  they  both  "good"?  (b)  If  (a)  is  true  which  is  better? 

The  second  set  of  questions  to  be  answered  has  to  do  with  speed. 

(a)  How  fast  are  the  routines?  (b)  Which  is  faster?  Finally,  is 
either  good  enough  and  fast  enough  to  be  highly  recommended? 

II.  Algorithms 
A.  RANDUL 

RANDUL  uses  the  algorithm  described  by  Lehmer*^  which  generates 
residues  of  some  Integer  by  a linear  transformation.  In  this  method 
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one  multiplies  the  current  random  integer  1 by  a constant  multiplier 
K and  keeps  the  remainder  after  division  by  m: 

new  I = K X old  I (mod  m) . 

This  is  termed  a linear  congruential,  or  more  specifically,  multi 
plicative  congruential  method  and  is  defined  formally  by  the  equation 


where 


and 


X ^ = ax  (mod  m) , n ^ 0 , 
n + 1 n ' ’ ’ 

x^  is  the  starting  integer  (seed),  x^  > 0 , 
a is  the  integer  multiplier,  a > 0 , 


m is  the  modulus,  m > x and  m > a . 

The  value  of  x^  should  be  odd,  and  the  value  of  a should  be  large  and 
odd.  Another  name  for  this  method  is  the  power  residue  method. 

By  choosing  m to  be  the  largest  integer  expressible  in  a single 
word  plus  one  (for  the  ASC,  2®^  + 1 = 21^7483649)  calculation  is  fast 
and  trivial.  The  residue  calculation  is  automatically  done  by  the  com 
puter's  overflow  handling  of  Integer  multiplication. 

Some  care  must  be  taken  in  choosing  a,  the  multiplier,  to  assure 
that  the  sequence  has  maximum  period  possible  for  a multiplicative 
generator.  The  maximum  period,  T,  for  m = 2®®  is  given  by® 

t(2'*')  = 2^  - 2 = = m/4  . 

Maximum  period,  although  desirable,  is  not  sufficient  for  randomness. 

A proof  given  by  Carmichael®  and  discussed  by  Knuth®  shows  that  the 


maximum  period  is  achieved  for  2”^  with  ^,  > 3 when  x^,  the  seed  is  odd, 
and  a = 3 or  5 (mod  8).  RANDUL  uses  the  multiplier,  2^®  + 3 = 65539 
which  meets  the  condition  65539  = 3 (mod  8). 

The  multiplier  65539  was  chosen  because  it  met  the  criterion  men- 
tioned above,  appeared  to  be  widely  used,  i.e. . found  in  much  document- 
ation,^ and  was  used  in  the  IBM-360' s RANDU.  Discussion  of  possible 
better  choices  will  be  covered  in  another  section. 

Once  the  power  residue  calculation  takes  place,  the  integer  is 
saved  as  a seed  for  the  (n  + l)-th  call  to  the  subroutine.  The  integer 
value  is  floated  and  multiplied  by  2“®^  to  distribute  the  values  on  the 
unit  interval. 

B.  VRANF 

It  is  convenient  to  consider  the  algorithm  of  VRANF®  to  be  com- 
posed of  two  sub-algorithms.  The  one  initializes  the  routine,  supplies 
the  first  1000  pseudo-random  numbers  that  are  requested,  and  provides 
the  starting  point  for  fulfilling  all  later  requests.  The  other  satis- 
fies all  later  requests.  Both  sub-algorithms  employ  integer  arithmetic, 
and  generate  odd  Integers  modulo  2®®  which  are  then  normalized  in  order 
to  supply  a finite  set  of  allegedly  uniformly  distributed  numbers  chosen 
from  the  open  interval  (0.0, 1.0). 

The  initialization  algorithm  is  based  upon  the  linear  congruential 
sequence  with  the  multiplier  1532i^550^7,  and  the  modulus  2®^.  The  se- 
quence does  not  meet  the  criterion  discussed  earlier  for  maximum  period, 
but  only  a small  sample,  the  first  1000  numbers,  is  used. 
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11.  From  an  Initial  positive  integer  which  is  supplied 
either  internally  or  by  the  user,  calculate  a seed  xi 
which  is  large,  positive,  and  odd. 

12.  Generate 

X , = ax  (mod  2®^) , 1 £ n ^ lOO  , 
n + 1 n 

and  remove  the  three  least  significant  bits. 

13.  X T=x  «-2x  + l,lSnS  100  , 

n,x  n n 

and  perturb  the  x^  in  order  to  ensure  that  all 
are  different  but  still  odd. 

14.  Generate 

X , = X,  ^ = a'^  ^x  , T (mod  2®^)  2 ^ k:S  10,  1 s n s loo  , 

and  remove  the  two  least  significant  bits. 

15.  Ensure  that  all  x . are  odd,  unique  and  computed 

11  ^ K 

modulo  2^®. 

At  this  point  1000  pseudo-random  integers  are  available  for  normal- 
ization and  transfer  to  the  output  array,  and  the  first  eight  entries 
are  stored  for  later  use  in  step  R3.  It  will  be  useful  to  denote  these 

^+1000'  1 ^ I ^ 8. 

All  remaining  requests  are  generated  according  to  the  following 
algorithm: 

Rl.  Choose  a pseudo-random  Integer  I 

I = 7 X 2"2®  ,1^15  8. 

R2.  X ♦-  X X (mod  2^®) , 1 s n s;  1000  . 
n n n + I 

^ + 1000  ^ 1 I ^ 8 • 

r4.  Normalize  x , 1 £ n ^ 1000  , 
and  transfer  to  output  array. 
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To  the  authors'  knowledge  there  is  no  rigorous  proof  that  the 
algorithm  given  by  R1-R3  produces  uniformly  distributed  odd  integers 
modulo  2^®.  Rather,  they  appeal  to  the  results  of  the  statistical 
tests  performed  on  the  normalized,  floating  point  results  obtained 
from  it. 


III.  Tests  of  Uniformity 

Statistical  tests  were  performed  on  RANDUL  and  VRANF  for  10  seeds, 
the  separate  default  seeds  and  nine  randomly  chosen  seeds.  A seed  must 
be  odd  and  no  greater  than  2147483647.  They  are  listed  in  Table  1,  and 
correspond  to  the  seeds  in  the  remaining  tables. 

TABLE  1.  Seeds  Used  in  Testing  VRANF  and  RANDUL 


Seed  No. 

Seed 

1 

1234 56781 

2 

379354235 

3 

246833 

4 

777654235 

5 

141597 

6 

73192881 

7 

429634011 

8 

2826305 

9 

1264396293 

10^ 

392396625 

lo’^ 

1009939513 

Default  seed  for  RANDUL. 
^Default  integer  for  VRANF. 


A.  Uniformity  of  singles  as  measured  by  the  X^-test 
, The  most  elementary  test  of  equidistributlon  is  to  partition 

the  unit  Interval  into  small  subintervals  of  equal  width  and  count  the 

I 

i 

number  of  terms  of  any  subsequences  that  fall  in  each  subinterval. ^ 

I 

The  values  are  then  compared  with  the  expected  values  by  means  of  a 
I chi-square  test. 

1 The  unit  Interval  was  divided  into  32  cells  using  sample  size 

1000,  64  cells  using  sample  size  5000,  and  128  cells  using  sample  size 
10000. 

Considering  Tables  2 and  3,  we  can  not  say  that  the  generated  num- 
bers are  not  uniform  by  singles  since 

var  (X®)  = 2 f , 

where  f is  the  nijmber  of  degrees  of  freedom,  and  in  any  one  column  at 
most  three  values  exceed  f + (2f)^/^,  which  is  within  acceptable  limits. 


[ 
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TABLE  2.  Uniformity  Test  Results  - RANDUL 


1000 

Nos.  32  groups 

5000  Nos.  64  groups 

lOOOO  Nos.  128  groups 

Seed 

100  P 

IBI: 

lOOPlX^SX^^) 

1 ^ 

1 ^ 

100P(X^fiX 

0 

■ 

30.7 

51.6 

64.4 

57.2 

53.5 

41.1 

89.4 

85.6 

96.9 

95.0 

B 

32.1 

58.9 

68.0 

68.8 

5.3 

B 

37.7 

81.0 

55.8 

27.2 

98.2 

2.7 

26.6 

30.9 

77.7 

90.0 

156.3 

96.0 

B 

41.0 

89.3 

70.0 

74.6 

129.8 

58.5 

7 

20.2 

6.8 

51.9 

16.1 

127.6 

53.1 

8 

32.0 

58.3 

56.6 

29.7 

159.7 

97.4 

30.8 

52.6 

71.7 

78.8 

127.5 

52.8 

32.3 

59.7 

59.5 

39.9 

125.7 

48.5 

TABLE  3.  Uniformity  Test  Results  - VRANF 


BHi 

1000 

Nos . 32  groups 

5000  Nos.  64  groups 

lOOOC  Nos.  128  groups 

Seed 

^0^ 

lOC  P (X2  s X 

0 

1CCP(X®SX 

0 

^0^ 

1 

20.4 

7.2 

63.0 

52.4 

133.  9 

68.0 

2 

24.6 

21.6 

58.8 

37.5 

98.2 

2.7 

3 

37.1 

79.0 

82.8 

95.2 

133.1 

66.2 

4 

21.6 

10.6 

76.1 

87.5 

132.6 

65.2 

5 

21.3 

9.7 

73.7 

83.2 

150.5 

92.4 

6 

41.3 

89.7 

65.6 

61.2 

129.5 

57.9 

7 

21.1 

9.1 

72. C 

79.6 

124.1 

44.2 

3 

34.8 

70.6 

68.2 

69.4 

109.5 

13.4 

9 

21.3 

9.7 

54.5 

23.1 

123.2 

42.1 

10 

18.2 

3.4 

52.0 

16.2 

136.2 

72.8 

I 
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B.  Iftiiformity  of  pairs  as  measured  by  the  x^'test 

Not  only  singles  should  be  uniform  on  a line  but  also  pairs 
should  be  uniform  over  a square.® 

Successive  pairs  of  pseudo-random  numbers  were  taken  as  points 
in  the  unit  square  which  was  divided  into  36  cells  for  5OO  pairs,  64 
cells  for  2500  pairs,  100  cells  for  5000  pairs.  A X®-test  was  then 
performed  on  each  of  the  distributions  obtained,  and  the  results  are 
given  in  Tables  4 and  5- 

The  probabilities  shown  in  the  third  and  fifth  columns  of  both 
tables  suggest  the  possibility  of  dangerous  geometrical  correlations 
of  successive  pairs.  On  the  other  hand,  the  last  column  in  each  table 
is  so  impressive  that  we  are  inclined  to  believe  that  both  generators 
are  satisfactory  with  respect  to  the  pair  distribution. 

Using  the  results  to  compare  VRANF  and  RANDUL,  we  can  not  say  that 
one  significantly  outperforms  the  other.  In  addition,  we  note  that  the 
results  with  the  default  seeds  are  at  least  as  good  if  not  better  than 
any  of  the  randomly  chosen  seeds. 
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TABLE  *)■.  Patrwlse  Uniformity  Test  Results  - RANDUL 


500  pairs  36  groups 

2500  pairs  64  groups 

5000  pairs  100  groups 

Seed 

Xo" 

100  P (X^ 

Hi 

lOOP(X^sXjj^) 

m 

100  P (X^  S X^2) 

1 

12.7 

69.6 

73.4 

98.5 

2 

86.9 

85.9 

97.1 

112.6 

3 

8.4 

51.2 

14.5 

94.7 

39.7 

k 

38.6 

68.8 

56.6 

29.9 

95.1 

40.7 

5 

33.2 

44.0 

83.7 

95.8 

91.9 

32.0 

6 

58.0 

99.1 

72.7 

81.2 

101.3 

58.2 

7 

30.1 

29.5 

61.0 

45.3 

123.7 

95.3 

8 

49.5 

94.7 

60.5 

43.3 

76.7 

4.7 

9 

35.2 

54.4 

51.8 

15.9 

109.4 

77.7 

10 

21.4 

3.5 

47.9 

7.9 

83.4 

13.1 

TABLE  5-  Pairwise  Uniformity  Test  Results  - VRANF 


500  pairs  36  groups 

2500  pairs  64  groups 

5000  pairs  loo  groups 

Seed 

100  P (X^  s X^®) 

Hi 

100  P (X^  S X^®) 

X ^ 
0 

100  P (X^  s X^^) 

61.0 

80.1 

92.8 

110.1 

79.1 

41.8 

59.4 

39.4 

121.1 

93.5 

24.6 

9.5 

71.1 

77.4 

100.2 

55.2 

37.1 

62.9 

90.1 

98.6 

86.3 

18.5 

33.2 

44.6 

47.6 

7.4 

81.1 

9.5 

38.6 

68.8 

63.3 

53.7 

94.4 

38.9 

41.1 

78.1 

65.7 

61.8 

124.2 

95.6 

8 

37.4 

64.1 

56.0 

27.8 

115.0 

87.1 

34.0 

48.1 

69.3 

72.6 

94.3 

38.5 

18.2 

0.9 

75.3 

86.1 

100.1 

54.9 
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C.  Uniformity  of  triplets  as  measured  by  the  X®-test 

Dividing  the  unit  cube  into  equal  cells  and  testing  for  uni- 
formly distributed  triplets  affords  another  test  of  equidlstribution. 
MacLaren  and  Marsaglia^  reported  a version  of  RANDUL  on  a 36  bit 
machine  with  the  multiplier  2^^  + 3 does  not  produce  a uniform  density 
in  three  dimensions.  Blood^°  explains  algebraically  how  the  triplets 
fall  on  a set  of  planes  in  the  unit  cube  predetermined  by  the  multi- 
plier and  the  modulus  and  with  RANDUL* s multiplier,  66339,  triplet 
correlations  will  always  be  found.  MacLaren  and  Marsaglia's^  test 
with  the  unit  cube  divided  into  1000  equal  cells  and  3 triples  taken 
from  each  set  of  10  numbers,  with  a sample  size  of  30,000,  was  performed 
on  VRANF.  The  results  are  given  in  Table  6. 


TABLE  6.  Triplet  Uniformity  Test  Results  - VRANF 


Seed 

X 2 
0 

loop  (X^^X^^) 

1 

998.9 

49.4 

2 

1012 

61.2 

3 

1112 

99.3 

4 

1009 

58.2 

5 

954.2 

15.3 

6 

1067 

93.2 

7 

1069 

93.6 

8 

1062 

91.6 

9 

1022 

69.4 

10 

964.9 

21.7 

■| 

■j 
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The  probability,  expressed  as  a percentage,  that  the  appropriate 
chi-square  variate  will  exceed  the  observed  test  value  may  be  occasion- 
ally large  or  small,  but  the  preponderance  of  large  values  leads  to 
the  conclusion  that  VRANF  also  falls  the  triplets  test. 

D.  Uniformity  as  measured  by  the  Kolmogorov- Smirnov 
goodness-of-f it  test 

Another  more  powerful  method  of  testing  equidistribution  is 

the  Kolmogorov- Smirnov  (KS)  goodness-of-fit  test.  The  chi-square 

tests  for  distribution  considers  observations  falling  in  a finite 

number  of  categories  while  the  KS  test  looks  at  the  random  real  number 

between  0 and  1 as  a continuous  distribution. 

The  statistics  K and  K are  calculated  by 
n n 

K^^  = /n  max  ’ 

- CO  < X < + OO 

and 

K = /n  max  (F(x)  - F (x))  , 

-oo<x<  + oo  ^ 

where  F (x)  is  the  observed  distribution  function  and  F(x)  is  the 
n 

expected  distribution  function.  Knuth^^  remarks  that  using  several 
tests  for  moderately  sized  n and  then  combining  the  observations  later 
in  another  KS  test  will  tend  to  detect  both  local  and  global  nonrandom 
behavior.  The  test  was  performed  with  groups  of  1000  observations 
1000  times  and  a final  KS  test  performed  on  these  results,  thus  using 
a total  sample  size  of  1,000,000. 


12 


TABLE  3.  Application  of  Kolmogorov- Smirnov 
Goodness-of-Fit  Test  - VRANF 
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If  we  choose  significance  levels  for  the  and  corresponding 
to  probabilities  greater  than  955^  or  less  than  and  count  the  fail- 
ures^ we  might  expect  a IC^  failure  rate.  Tables  7 and  8 show  the 
results  of  the  1000  KS  tests  with  local  failures  falling  within  the 
expected  values.  However,  the  result  of  the  KS  test  applied  again  to 
the  1000  individual  test  results  detects  a global  nonrandomness  using 
seeds  6 and  7 in  VRANF  and  seed  7 in  RANDUL.  The  KS  test  evinces 
RANDUL  meets  the  criterion  for  randomness  better  than  VRANF. 


E.  Runs  test 

A strong  test  recommended  by  Knuth,^^  strength  being  implied 
by  the  fact  that  a large  number  of  random  number  generators  fail  it,  is 
a test  for  "runs  up"  and/or  "runs  down".  As  implemented,  the  test  tab- 
ulates "runs  up"  counting  six  categories  of  runs  of  lengths  1 through 
5,  and  runs  of  length  6 or  greater.  Since  adjacent  runs  are  not  inde- 
pendent the  results  are  evaluated  by  the  statistic 


V = — I (COUNT(i)  - nb.)  (COUNT(j)  - nb,)a^. 


where  a^j  and  b^  are  tabular  covariance  values  derived  and  calculated 
in  Ref.  (12).  The  statistic  V then  has  a chi-squared  distribution  with 
6 degrees  of  freedom.  The  results  are  presented  in  Tables  9 and  10. 

RANDUL  fails  the  runs  test.  Of  the  50  entries  in  Table  9,  sixteen 
exceed  0.9,  which  indicates  a highly  significant  deviation  from  the 
hypothesis  that  the  numbers  are  uniformly  distributed.  VRANF  survives 
much  better  with  only  8 of  50  entries  exceeding  0.9  which  we  might  term 
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TABLE  9.  Probabilities  that  is  Less  Than  the  Observed 


Value  for  RANDUL 

Subjected 

to  the  Runs 

! Test 

Gpl 

Gp2 

Gp3 

Gp4 

Gp5 

0.927 

0.979 

0.903 

0.771 

0.233 

0.556 

0.217 

0.875 

0.892 

0.154 

0.296 

0.482 

0.799 

1.0 

0.532 

0.027 

0.091 

0.998 

0.190 

0.652 

0.825 

1.0 

0.267 

0.324 

0.203 

0.826 

0.996 

0.907 

0.303 

0.722 

0.781 

0.743 

0.782 

0.820 

1.0 

0.996 

0.409 

0.875 

1.0 

0.967 

0.569 

0.995 

0.780 

0.993 

0. 925 

0.095 

0.115 

0.095 

0.916 

0.277 

*The  default  sequence  for  which  the  group  size  is  4133 
numbers.  For  all  other  seeds  the  group  size  is  5133 
numbers . 


TABLE  10.  Probabilities  that  is  Less  Than  the  Observed 
Value  for  VRANF  Subjected  to  the  Runs  Test 


Gpl 

Gp2 

Gp3 

0.698 

0.034 

0.969 

0.355 

0.536 

0.083 

0.216 

0.458 

0.260 

0.557 

0.773 

0.418 

0.352 

0.565 

0.505 

0.622 

0.728 

0.942 

0.941 

0.349 

0.132 

0.444 

0.946 

0.701 

0.664 

1 

0.457 

0.800 

0.862 

0.226 

0.998 

Gp4 


0.071 
0.489 
111 
985 
154 


Gp5 


0.347 

0.917 

0.504 

0.120 

0.489 

0.761 

0.296 

0.981 

0.566 

0.588 


*The  default  sequence  for  which  the  group  size  is  4133 
numbers.  For  all  other  seeds  the  group  size  is  5^33 
numbers. 
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as  not  showing  a significant  deviation  from  the  hypothesis  that  the 
numbers  are  uniformly  distributed. 

F.  Difference  test 

In  the  absence  of  a rigorous  proof,  the  author  of  VRANF  elicits 
its  validity  by  its  performance  on  a test  based  upon  differences. ® 

If  we  have  two  statistically  independent  random  variables  distri- 
buted uniformly,  xi  and  Xg,  the  difference  variable  yi  = xg  - x^  has  a 
triangular  distribution,  and  the  variable 

2i=  yi  + 1 yi  ^ 0 

= yi  yi  > 0 

again  has  a uniform  distribution.  Using  this  principle,  we  may  test 
the  outputs  of  a random  number  generator,  x^,  x^,  x^^,  x^^^^  = x^, 

their  first  differences,  Zi,  Zg,  •••,  = zi,  the  first  differ- 

ences of  the  z-variables,  etc,  for  uniformity. 

A sample  of  10000  numbers  was  chosen  to  start,  the  first  50  differ- 
ences of  the  x-variables  were  calculated,  and  uniformity  was  tested  by 
means  of  the  X^-test  with  the  unit  interval  divided  into  25  cells. 

Figure  1 shows  plots  of  the  theoretical  X^-dlstribution  for  2h 
degrees  of  freedom,  and  the  observed  results  obtained  from  the  default 
sequences  of  RANDUL  and  VRANF.  It  Illustrates  that  VRANF  remains  stable 
for  50  differences,  and  the  break  in  the  curve  for  RANDUL  at  a percentile 
value  of  about  kO  illustrates  that  its  randomness  collapses  at  the 
twenty-second  difference.  VRANF  consistently,  for  all  seeds,  matched  the 
expected  chi-square  distribution  for  2h  degrees  of  freedom  while  RANDUL 
failed  at  the  twenty-second  difference  for  all  seeds. 
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Fig.  1 — Chi-square  distribution  of  differencing  test  for  defauit  seeds  for  24  degrees 
of  freedom.  Results  are  shown  for  theoretical  distribution  and  the  observed  results 
obtained  from  RANDUL  and  VRANF. 


G.  Spectral  test 

Coveyou  and  MacPherson^  formulated  a test  to  measure  the  sta- 
tistical Independence  of  successive  n-tuples  of  random  numbers  using 
finite  Fourier  analysis.  In  a proof  detailed  In  their  paper,  It  Is 
shown  that  the  Fourier  coefficients  are  sufficient  statistics  them- 
selves, and  a computational  method  Is  presented  for  calculating  the 
coefficients  over  the  whole  period  for  a linear  congruentlal  generator. 
By  comparing  the  Fourier  coefficients  of  the  transform  of  a truly  uni- 
form random  sequence  with  those  of  the  transform  of  a uniform  random 
number  generator  for  given  n-tuple  sizes,  the  spectrum  of  deviations 
of  the  generator  output  from  uniform  randomness  Is  obtained.  For 
linear  congruentlal  sequences,  the  Fourier  transform  over  the  whole 
period  has  the  form 


f(si,  ...,s  ) ■ — L exp 
” “ OSk<m 


= — S exp 
0^k<m 


where 


= s^  + Sg  a + 83 


_n  - 1 


and  a Is  the  multiplier.  Assuming  the  sequence  has  maximum  period  this 
reduces  to  the  sum  of  a geometric  series 

where  0(x)  a 1 If  x Is  an  Integer,  and  0 otherwise.  The  quantity 
f(si,  therefore,  either  vanishes  or  has  unit  magnitude. 
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We  may  physically  Interpret  f(si,  s^) /m  as  the  amplitude  of  the 
n-dlmenslonal  plane  wave 

t ) = exp  2n  1(4^  ti  + ...  + — t)  , 
n I m m n I 

and  assign  a wave  number  t>  corresponding  to  Its  frequency 


V = / si^  + . . . + 8 ^ when  I s ^ — 

n ' k'  2 • 

For  linear  congruent lal  sequences  of  maximum  period,  the  smallest 
nonzero  wave  number  In  the  spectrum  Is  given  by 


= min /si'^  + > 


where  the  minimum  Is  taken  over  all  n- tuples  of  Integers  excluding 
(sx,  S2,  . . . , s^)  = (0,  0,  . . . , 0)  subject  to  the  constraint  s(a)  = 0 (mod  m) , 
To  rate  the  effectiveness  of  a multiplier  In  a linear  congruentlal 
sequence  of  maximum  period  we  calculate 

n/2 

c = I . 

" (n/2)  .'m 

Values  of  of  order  unity  or  larger  correspond  to  randomness, 
small  values  correspond  to  nonrandomness.  Large  values  Imply  large 
wave  number  (short  wavelength)  deviations  that  will  tend  to  average 
out  over  the  entire  cycle.  Long  wavelengths  Imply  a systematic  de- 
parture from  uniformity. 

Implementing  the  spectral  test  Involved  using  the  ASC  extended 
precision  package  for  integer  arithmetic,  SAC-1,  since  the  computer 
algorithm  described  in  Ref.  (13)  requires  exact  Integer  operations. 
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The  author  of  Ref.  (13)  qualifies  the  algorithm  by  stating,  "the  method 

will  get  into  an  infinite  loop  in  some  cases."  Modifying  the  algorithm 

to  avoid  infinite  loops  leads  to  a very  long  search  for  the  quantity 

. Since  the  extended  precision  package  is  implemented  by  means  of  a 
n 

linked  list,  it  is  very  space  and  time  consuming,  and  a lengthy  search 
may  be  prohibitive.  This  was  the  case  in  solving  for  coefficients  C3, 
C4,  and  C5  for  VRANF  as  indicated  in  Table  11. 


TABLE  11.  Spectral  Test  Results 


Coefficients 

RANDUL 

VRANF 

Cl 

2.0 

2.0 

C2 

3.1412 

2.19 

C3 

0.0000025 

•it- 

C4 

0.00003092 

* 

C5 

0.000355 

* 

*Non- converging  algorithm. 


From  results  shown  in  Table  11  RANDUL  must  be  rejected  as  a gener- 
ator where  the  uniformity  of  distribution  of  n-tuples  greater  than  2 
is  a critical  factor. 

The  statistics  shown  for  VRANF  are  not  as  potent  since  the  spec- 
tral test  rates,  over  the  full  period,  the  multiplier  of  a linear- 
congruentially  generated  sequence,  and  the  congruential  algorithm  is 
only  used  to  generate  the  first  1000  numbers  in  VRANF.  A true  rating 
of  VRANF  by  a theoretical  test  is  non-trivial  since  the  full  period 
has  not  been  mathematically  predicted. 
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F.  A.  Blood^°  discusses  and  predicts  the  harmful  correlations 
found  by  the  spectral  test  in  the  power  residues  generated  by  RANDUL. 
He  shows  algebraically  how  to  pick  a multiplier  to  eliminate  the 
problem  of  geometrical  correlations  of  n-tuples  3 or  greater  and 
supplies  a set  of  multipliers  which  meet  his  criterion. 


H.  Autocorrelation  test 

A quantity  of  interest  prior  to  employing  pseudo-random 
numbers  ir  many  applications  is  the  autocorrelation  function  of  such 
a sequence.  Several  definitions  of  this  quantity  exist,  and  for  the 
present  purpose  it  is  convenient  to  choose  one  that  bounds  it  by  unit 
magnitude.  The  serial  correlation  coefficients  provide  such  a defi- 
nition and  are  given  by 


where 


CU) 


E[(Xn  - - X)] 

E[(x^  - x)2] 


N 

E[Z^]  = (1/N)  E Z^, 
n = 1 


and  X = E[x^],  with  the  x^  being  the  pseudo-random  numbers  that  have 
been  generated  by  the  routine  under  test.  The  convolution  appearing 
in  the  numerator  has  been  calculated  by  first  obtaining  its  finite 
Fourier  transform  with  the  routine  FOURTH, which  reduces  the  con- 
volution to  multiplication,  and  then  inverting  the  result.  With  this 
technique,  the  values  of  C(-t)  have  been  obtained  from  the  first  2^“^ 
numbers  generated  by  RANDUL  and  VRANF  for  the  seeds  given  in  Table  1, 


and  for  0 S -t  ^ 2^^.  The  sample  means  and  the  values  of  C(^)  for 
t = 1,  2,  4,  8,  l6  are  given  in  Tables  12  and  13  for  RANDUL  and 
VRANF  respectively. 

Various  assumptions,  all  of  which  involve  the  normal  distri- 
bution at  some  stage,  about  the  quantities  for  which  the  serial  cor- 
relations are  of  interest,  indicate  that  for  t > 0 and  large  n,  the 
C(^)  tend  to  be  normally  distributed  with  mean  and  variance  0(1/ n) . 
With  the  assumption  that  this  theory  is  applicable  to  the  serial  cor- 
relation coefficients  of  a sequence  obtained  from  the  uniform  distri- 
bution, the  probability  of  obtaining  jcC-Oj  greater  than  1.6  n"^/^  = 
0.0125  is  10^,  and  that  of  obtaining  a result  greater  than  2n"^^^  = 

0.0156  is  55^. 

Of  the  fifty  entries  in  the  last  five  columns  of  Table  12,  there 
are  three  greater  than  the  5^  critical  point,  and  four  greater  than 
the  Ic^  critical  point.  Similarly  in  Table  13,  the  corresponding 
numbers  are  two  and  four.  These  results  indicate  that  both  ElANDUL 
and  VRANF  have  acceptable  serial  correlation  properties. 
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TABLE  12.  Results  of  the  Autocorrelation  Test 
Applied  to  RANDUL  and  Sample  Mean 
Values 


Seed  No. 

X 

C(1) 

C(2) 

C(U) 

— 

C(8) 

C(i6) 

1 

C.^93l 

-0.0102 

0.0020 

0.0099 

0.0170 

0.0200 

C.5Cl2 

-0.0009 

0. 00^5 

0. 0051 

0.0001 

-0.0066 

5 

C.5C27 

-0.0035 

-0.0051 

0.0052 

-0.0159 

0.0060 

C. 5009 

0.0055 

-0.0121 

-0.0057 

-0.0136 

0.0086 

0. 5005 

-0. 0057 

0.0085 

l.0006 

-0.0080 

-0.0100 

r 

c 

0. 5055 

-0.0027 

-0.0017 

0.0022 

0.0018 

0.0018 

1 

C.^93l 

-c.coai 

C.QQl'k 

-0.0102 

-0.0057 

-0. 0060 

3 

0.5005 

-0.0050 

-0.0107 

-0.0088 

0,0020 

-0.0121 

o 

0.5026 

0.0012 

-0.00^5 

-0.0104 

0.0014 

0.0000 

1C 

O.i^-997 

0.0116 

0.0019 

-0.0089 

-0.0024 

0.0007 

TABLE  15.  Results  of  the  Autocorrelation  Test 
Applied  to  VRANF  and  Sample  Mean 
Values 


Seed  No. 

X 

C(l) 

C(2) 

C(4) 

C(8) 

1 

c(l6) 

_L 

0.4952 

-0.0027 

-0.0207 

-0.0009 

0.0105 

0.0069 

2 

0.499c 

-C.CC4o 

0.0000 

0.0069 

-0.0184 

-0.00I2 

3 

0. 5006 

0.0005 

-0.0036 

0.0115 

0.0085 

0.0117 

4 

0.5003 

0.0056 

-0,0070 

-0.0092 

-0.0086 

-0.0032 

i; 

y 

0.5C44 

0.0015 

-0.0113 

-0.0029 

-0.0090 

-U.C030 

6 

0.4932 

-0.0157 

C.CC52 

-0.0129 

0.0083 

-0. 0020 

7 

0.5011 

-0.0027 

0.0092 

-0.0001 

-0.0012 

-0.0056 

3 

0.4986 

-0.0028 

-0.0070 

-0.0124 

-0.0120 

0.0115 

9 

0.4991 

0.0079 

0.0001 

-0.0005 

-0.0002 

0.0005 

10 

0.5019 

0.0045 

-C.014C 

-O.OlOO 

0.0049 

0.0163 

IV.  Timing  Tests 


The  subroutine  RANDUL  must  be  called  once  for  every  number  gener- 
ated. The  routine  VRANF,  however,  uses  the  machine  vector  capability 
to  generate  and  store  an  array  of  1000  numbers  regardless  of  the  value 
N,  the  number  of  numbers  requested.  The  N numbers  are  then  moved  to 
the  output  array.  Subsequent  calls  generate  a new  array  of  1000 
numbers  when  the  original  1000  are  depleted.  Thus,  time-wise  use  of 
VRANF  is  optimized  when  N = 1000. 

Results  of  tests  in  Table  14  show  VRANF  is  faster  than  RANDUL  by 
a factor  of  1000  when  the  former  is  used  at  its  maximum.  However, 
RANDUL  is  a more  efficient  generator  when  only  a few  numbers  are 
desired.  Timing  was  performed  by  routine  SECOND. 


TABLE  14.  Timing  Tests  RANDUL  vs.  VRANF 


Test 

No. 

Test  Description 

RANDUL 
time  in  sec. 

VRANF 

time  in  sec. 

1 

Default  seed  - Single  call  N = 1 

0.31  X 10“^ 

0.95  X 10-3 

2 

Seed  initialization  only 

0.26  X lO""^ 

0.93  X 10-3 

3 

Average  time  - 1000  calls  N = 1 

0.32  X 10"^ 

0.44  X 10' 

4 

N = 1000® 

0.20  X 10'^ 

0.21  X 10'^ 

5 

VRANF  - Single  call  N = 1001 

0.41  X 10'^ 

6 

VRANF  - Single  call  N = 3000 

0.63  X 10" 

a 

Includes  program  looping  for  RANDUL. 
Single  call  to  VRANF  with  N = 1000. 


■.—I  • • 

- 


24 


V.  ReconsnendaClons  and  Restrictions 

We  conclude  that  VRANF  outperforms  RANDUL  both  in  speedy  period, 
and  somewhat  In  statistical  performance,  although  both  have  difficulty 
uniformly  filling  the  unit  cube. 

RANDUL  is  recommended  for  usage  when:  1.  generating  a small 

mmiber  of  numbers;  2.  the  user  desires  the  property  that  an  individual 
number  will  not  repeat  until  the  whole  sequence  repeats;  3.  the  user 
desires  random  integers  of  range  1 to  2^^  + 1. 

VRANF  is  recommended  for  general  usage,  however,  the  user  might 
note  a few  potential  problems  discovered  in  using  it.  The  seed  value 
supplied  in  the  subroutine  call  is  altered  by  VRANF  and  the  altered 
value  returned.  Calling  VRANF  without  a call  to  SETVR  for  initializing 
a seed  presently  does  not  generate  the  same  sequence  as  calling  SETVR(O) . 
VRANF  documentation,  however,  purports  this  to  be  true.  The  routine 
will  be  changed  to  generate  duplicate  sequences  in  this  case,  in  order 
to  conform  with  existing  documentation.  To  use  the  subroutine  effi- 
ciently the  user  is  responsible  for  dispersing  the  supplied  array  of 
random  numbers,  and  requesting  more  when  the  queue  is  empty,  which 
may  involve  changing  the  logic  of  existing  routines  currently  using  a 
random  number. 

Neither  routine  can  be  used  for  producing  matching  sequences  on 
the  ASC  and  CDC.  There  exists  a slow  random  number  generator,  RANF, 

for  this  purpose  which  is  recommended  for  use  only  for  the  duplicate 
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sequence  problem. 
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Appendix  A.  RANDUL  - Availability  and  Source  Listing 


A.  RANDUL  resides  on  the  SPL  whose  objects  are  currently  on  SYS.OLIB. 


SUBROUTINE  RftNOUL ( lY tVFL ) 

CD  1.  IDENTIFICATION  NARE:  2.  CL ASS  IF IC AT ICNT  3.  DATE: 

CO  RANDUL  C5  93C^76 

CO  4.  title: 

CO  UNIFORMLY  DISTRIBUTED  RANDOM  NUMBER  OENTRATCP 

CD  S.  AUTHORCS):  ORG AN  I Z AT  ION (S ) : 

CO  UNKNONN  : OAKRIOGE  LABORATORY 

CO  MODIFIERCS):  CRCANIZATIONCS): 

CO  TERRY  CARUTMERS;  NRL  RCCT  G.H.  HERLINGJ  NRL  STD 

CO  CONVERTER(S);  ORG ANIZ AT ION(S ): 

CO  TERRY  CARUTMERS;  NRL  RCC 

CO  reviewer:  organization: 

CO  SIDNEY  PRAHL;  TEXAS  INSTRUMENTS.  INC. 

CO  6.  DESCRIPTION/PURPCSE  (SPECIFY  SUBROUTINE,  FUNCTION,  PR  PROGRAM): 

CO  SUBROUTINE  GENERATES  A UNIFORMLY  OlSTRieUTEO  PScUOC-R ANOOM 

CD  NUMBER  ON  THE  INTERVAL (0 .0 , 1 .0 ),  END  POINTS  EXCLUDEO,  AND/Ok 

CO  INTEGER  ON  THE  INTERV At <0 , 2** 3 1- 1 ),  USING  LINEAR  CeNGRU=NTIAL 

CO  ALGORITHM.  THE  USER  MAY  SUPPLY  A STARTING  SEED  BY  A SINGLE  CALL 

CO  TO  RANOUS(ISEEO)  BEFORE  CALLING  RANDUL,  CR  CALL  RANDUL  ALCNt 

CO  WHICH  SUPPLIES  A DEFAULT  SEED.  SUCCESSIVE  CALLS  TO  RANJUL  CAN 

CO  GENERATE  A NONREPEATING  SEQUENCE  OF  MAXIMUM  LENGTH,  2**29. 

CO  BASIC  USE  IS  FOR  QUICK  GENERATION  OF  SHORT  SEQUENCES  OF  PSEUDO- 

CO  RANDOM  NUMBERS  THAT  APPEAR  TO  HAVE  BEEN  RANDOMLY  SELECTED  ERO"* 

CD  THE  UNIFORM  PROBABILITY  DENSITY. 

CD  7.  CALLING  SEQUENCE  OR  OPERATIONAL  PROCEDURE: 

CO  CALL  RANOULdY.YFL) 

CO  8.  ARGUMENTS  (TYPE  AND  SIGNIF ICANCE ) AND/Uh  INITIAL  CONDITIONS: 

CO  lY- INTEGER, OUTPUT  RANDOM  INTEGFD 

CD  YFL-FLOATING  PT,  OUTPUT  RANDOM  NUMBER 

CO  ISEEO-INPUT  INTEGER,  I TO  10  DIGIT  ODD  INTEGER  .LE . 

CD  (2147483647)  WHICH  DETERMINES  THE  UNIQUE  SEtUINCE. 

CO  NOTE:YFL=NORMALIZcD  IY. 

CO  9.  ENTRY  POINTS  (WITH  ARGUMENTS): 

CO  CALL  RANDUL(IY,YFL) 

CD  CALL  RANOUSdSEED) 

CO  10.  LANGUAGE/LINITATIONS/ACCURACY/SYSTEM  DEPENDENCES: 

CO  RANDUL  IS  A MULTIPLE  ENTRY  SINGLE  PRECISION  FORTRAN  IV  SU3R0UT1K: 

CO  AND  MUST  BE  COMPILED  UNDER  THE  C OPTION  OF  THE  sx  COMPILER  AS  iHt 

CO  ALGORITHM  DEPENDS  UPON  OVERFLOW  AND  UNDERFLOW  TO  SUPVIVc.  HIS 

CO  SUBROUTINE  IS  TRANSFERABLE  ONLY  TO  ANOTHER  32  BIT  MACHINE  WHICH 

CO  ALLOWS  OVERELCW/UNOERFLOW  AND  RETAINS  THE  LEAST  SIGNIFICANT  PART 

CD  OF  THE  INTEGER  MULTIPLY  OVERFLOW.  THE  SUBROUTINE  ALSO  NEEDS  A 

CO  CONSISTENT  ENTRY  POINT  STRUCTURE. 

cu  II.  routines  called: 

CO  NONE 

DATA  IX/392396624/ 

IY=IX*65939 
IF(IY)  9,6,6 

5 IY=IY42147483647*1 

6 YFL^IY 
IX  = IY 

YFL  = YFL*  .469661 3E-3 
RETURN 

ENTRY  RANUUS(IX) 

RETURN 

END 
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Appendix  B.  VRANF  - Availability  and  Source  Listing 


A.  VRANF  resides  on  SYS.OLIB. 

bUBROUTlNe  VRANF  (RA«LFN) 

CO  1.  ICtNT  IF  ItAI  ION  NANE:  ClASS  IF  ICAT  I -N:  ?.  oAT'': 

CO  VRANF 

CO  4.  iiile: 

CO  UNIFORNLT  OISIRIBUTEC  VECISR  kA.NOCF  NUFB'R  G'^N^RMSK 

CO  b.  AUIHSKCSi:  9RGANI/AIMN(S): 

CO  OICA  KOIH-II  AUSTIN 

CO  XOOIFIERCS):  ORGAN  WATIONCS>J 

CO  NONE 

CO  CONVERTERtS):  OKG AN i Z AT 10N< S ) : 

CO  lEKKY  CARUTHFRS  RCC-ARL  G.H.  « t Rl T NG t R TC - NRL 

CO  RhVibNEK:  ORGANIZATION: 

CO  SIDNEY  PHAHL:  TEXAS  TNSTRUNE NT S . INC. 

CD  6.  OESCRIPinN/PURPOSE  ESPECIpy  SUi'ROUT  INE  . FUNniTN,  OR  Rk'GRA*); 

CD  IHIS  ROUHNt  USES  VECTOR  SOEEO  TO  GcNERATE  AN  ARRAY  Or  RANOO"' 

CO  NUNBFRS  UNIFORNLY  OISTHIBUTEC  ON  T><E  RANGE  ( i . f . 1 . ' ) , l Nt  POINTS 

CO  excluded,  an  INIHiCATE  ALGORITMti  C^'RAINING  LINEAR  G ONC-RUENT I Al 

CO  GENcKAIICN  NIIM  RANDOM  CROSS  MULTIPLICATION  IS  USED.  US^B  MAY 

CD  SUPPLY  A STARTING  SEED  BY  A SINGLE  CALL  TO  j>'TVR<INIT>  oEEORE 

CO  calling  VRaNF, otherwise  a RUILT  IN  value  iS  USED.  ArTt»  ANY  Call 

CO  VRANF  MAINTAINS  liO:  PARALLEL  STREAMS  fF  RANDOM  NUMBERS  AND 

CO  REGtNERAIES  ONLY  WM£N  THES'  ARE  USED 

CO  T.  CALLING  SE'JUENCc  OR  OPERATIONAL  PROCECIkE: 

CO  call  VR ANF(RA,LEN)-STAN0AR0  SUBROUTINE  ENTRY 

CO  CALL  S£ IVHtINITJ-tNTRY  POINT  FOR  INITIAL  SEED  MOD  IF TC AT  I ON. 

CO  8.  ARGUMENTS  tTYPE  AND  SIGNIFICANCE)  AND/CH  INITIAL  CONDITIONS: 

CO  RA-REAL-STARTING  ADDRESS  OF  OUTPUT  RANCOP  ARRAY 

CD  LtN- INTEGER-LENGTH  OF  ARRAY  TO  BE  fillED 

CO  INII-TNIEGER-  I TO  lii  DIGIT  ODD  INTEGER  .LE.  2*  • ? I- U E I ••  TO  E 360  7i 

CO  STARTING  SEED  WHICH  0£TtRMIN'’S  UNIQU=  SttC^NCE. 

CO  V.  ENTRY  POINTS  (WITH  ARGUMENTS): 

CO  CALL  SEIVKLINID-USED  AS  A ONE  TIME  ONLY  CALL  TO  GEN^kAIE  A 

CO  UNIOUr  SEQUENCE  OF  RANDCM  NUMBERS.  ZERO  AS  A 

CO  SEED  GENERATES  THE  OEFAILT  SERIES.  WITm  NO  CALL 

CO  TO  SETVR  VRANF  SUPPLIES  A SE  EO.  1 08  2 1 Z 300  7 f mHICa- 

CO  ALWAYS  PRODUCES  Th£  SAME  SEQUENCE  OF  P£RIOD 

CO  ?«*27000. 

CO  ARGUMENT-INIT  -1-lC  DIGIT  ODD  INTEGER 

CO  10.  LANGUAGE/LIMTIATIONS/ACCURACY/SVSIEM  DEPENDENCES: 

CD  ROUIINE  IS  A MULTIPLE  ENTRY  SINGLE  PRECISION  FORTRAN  IV  SUB- 

CO  ROUTINE  ANO  MUST  BE  COMPILED  UNDER  C OPTION  OF  NX  COMPILER,  AS 

CO  THE:  ALGORITHM  DEPENDS  UPON  OVERFLOW  ANO  UNDERFLOW  TO  SURVIVc. 

CO  COMPILING  UNDER  X OPTION  PRODUCES  A VEflCR  HAZZARO  T“AT  Mfty 

CO  safely  ignored  under  l OPTION.  ROUTINE  IS  TRANSEERAdLE  ONLY  To 

CO  ANOTHER  32  BIT  MACHINE  WHICH  ALLOWS  mvER'LOW/UNOERFlOw  aND 

CO  IN1EGER  MULITPLY  OVERFLOW.  ALSO  NEEDS  CONSISTENT  ENTRY  POINT 

CD  STRUCTURE.  ROUTINE  IS  WRITTEN  WITH  II  VECTOR  PHILOSOPHY  IN  MIND. 

CO  ADVANTAGES  OF  SPEED  HILL  BE  LOST  IF  EXECLTEO  AS  SCALAR  CODE 

CD  11.  ROUTINES  CALLED: 

CD  NONb 

DIMENSION  JR(10D8),JRP(100,10),RA(LEN) 

EQUIVALENCE  (JR,JRE) 

DATA  ISTARI,  ISEEO,  JANO,  KNOR ■ , JSE T , JR E M , TN9R m 
• / Z3C327A39,Z4BS768e7,Z«FFFFFFF,Z3AlOO?0:,  CtlCOO,  u.J  > 

C 

C RETURN  WITH  NO  ACTION  IF  LEN  < 1. 

IF  ClEN  .LT.  1)  RETURN 
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oo  no  on  n n n n n n n on  on  nr  nn 


INIIIALIZE  FIRST  lOOC  NOS.  THIS  HAS  NOT  BEEN  DONE. 

IF  CJSET  .EU.H)  GO  TO  n 
8C  NEED  = LEN 
90  HEAD  - NEED 
KREM  = JREM 


00  1 HAVE  ENOOr.N  UNUScO  NUMn£>(S  TO  hiLL  THIS  Rcu'JtST? 

IF  AjKEH  .Gi.  need;  GO  TO  n,; 

NO.  BUI  I'LL  USE  mHAT  T HAVE  AVATt^bLE  FCH  Nv*«. 

NUN  = JRE> 

NEEP  = NEED  - JRE" 

JREN  - C 

GO  10  nc 

VES*  1 HAVE  ENOUGH. 

100  NUN  * NEED 

JREN  = JKEN  - held 

NEED  = 0 

FLOAT  NUNBEKS  AND  NOVE  TO  OUTPUT  ARhAt. 

119  00  120  H = liNUN 

129  R ALN*LEN-N£AO>  = JRtN*l'iC  i-RREN ) 

010  1 NOVE  ALL  THAT  WERE  KEwUESTEu? 

IF  LNEEO  .GT.  0)  GO  T"  l«»l 

YES.  REQUEST  IS  FILLEP.  NOW  NORNALU'  TO  RANGE  ('.C.l.-)  6V 
NULTTPLICATION  by  RNOKN. 

00  130  N = ItLcN 

130  RALNl  = KNORNARALNl 
RETURN 

NO.  1 OlON'l  have  enough.  $0  N^h  i GENERATE  NEXT  lOf'  No"uckS. 
140  JP  s JRLOTTl* TNORN  « 1 
DO  ISO  J > UlOOG 

ISO  JR(J)  > and  (JR(J)«JR( J4JP). JANO) 

00  IbO  J = 1*8 
160  JR(J«100('l  = JR(J) 

JREN  > U09 

NON  GO  BACK  SNO  CONPLCTE  REST  OF  REQUEST. 

GO  10  90 

JUNP  TO  HERE  IF  VRNON  CALLEO  WITHOUT  PRIOR  CALL  TO  SETVR. 

10  IFLG  = 1 

INT  3 ISTART 
GO  TO  20 
C 

C ••••«••••« 

c 


30 


o r> 


I 

i 


ENTKV  SEIVK  UNIT) 

INT  = (tNlT«7b6R3I)*UNir«9343(i7^)  * ISTAPT 
IFLO  = 0 
?0  JSEl  = I 
JRE»  = ICO? 

TNOKM  = 7*KNflRM 
JRllJ  = 2*INr*ISlAhT  ♦ 1 
C 

C GENEKAIE  HkST  1C3  KAhOOH  NU»<oEkS  ST.ALAH  LOOP. 

J = 2 

30  J«LJ>  = LSHFl TSElD* JRCJ-1 )«-3i 
J = J ♦ 1 

IF  (J.LT.131)  GO  TO  3'^i 
00  32  J = UlSf 
32  JRLJ3  = 2*JRlJ>  ♦ 1 
C 

C ROOIFY  NUHBERS  TO  AVOID  RULITPLE  GUPLICATt  3't.UENCEl) 
DO  3!>  J = 1,100 
35  JRIJO  -■  J«LJ>  ♦ 15<»fl63t)2«J 

USE  VEGTORT/EO  loops  to  extend  seed  10  t.  = NGTF  Of 
00  55  K = 2,i: 

IF  LJKEM.LI.OO  CO  TO  55 
DO  50  J = 1,1P» 

50  JREtJ.K)  = LSMF(iStEU*JRE(J»K-l),-2> 

55  CONIINUE 

CO  57  J = l.lOilo 
57  JR(J)  = OKtJRU),!) 

C ROOIFY  NURStRS  TO  AVOID  RULTTPLE  DUPLICATE  SFiJUCNCS 
DO  50  J = Iflwj;' 

60  JR(J)  = JPLJ)  ♦ 176S535*J 
00  73  J = 1,133: 

70  JRLJ7  = ANDlJPLJi,JAftn  ) 

DO  72  J = l,tJ 
72  JR(J»nOO)  = JR(J) 

C 

IF  LIFLG.EJ.J)  MEMIF^ 

GO  10  80 
C 

END 


